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We report on combined neutron and resonant x-ray scattering results, identifying the nature of the 
spin-orbital ground state and magnetic excitations in LuVOa as driven by the orbital parameter. In 
particular, we distinguish between models based on orbital Peierls dimerization, taken as a signature 
of quantum effects in orbitals, and Jahn-Teller distortions, in favor of the latter. In order to solve 
this long-standing puzzle, polarized neutron beams were employed as a prerequisite in order to solve 
details of the magnetic structure, which allowed quantitative intensity-analysis of extended magnetic 
excitation data sets. The results of this detailed study enabled us to draw definite conclusions about 
classical vs quantum behavior of orbitals in this system and to discard the previous claims about 
quantum effects dominating the orbital physics of LuVOa and similar systems. 

PACS numbers: 75.70.Tj; 75.10.Dg; 78.70.Nx; 78.70.Ck; 75.40.Gb 


Geometrically frustrated magnetism with its crucial 
role of quantum effects in purely spin-systems is a 
well-established held with new, exotic phases emerging 
[1]. Gan we, however, have a similar situation in a 
completely different context, namely, in orbital physics? 
A positive answer to this question would open up a new 
held and dehne a new class of materials. 

The possible role of quantum huctuations in orbital 
physics is a very interesting and important question. For 
small objects such as Jahn-Teller (JT)-active molecules 
or isolated JT impurities in solids, such quantum effects 
are well known and constitute a big held of vibronic 
effects in JT physics [2, 3]. On the other hand, in 
concentrated solids we practically always ignore these 
effects and treat orbitals (quasi) classically. Therefore 
all the more exciting were the suggestions [4, 5] that 
orbitals may behave as essentially quantum objects, in 
particular in some perovskite vanadates [ 6 , 7], up to the 
formation of orbital singlets in YVO 3 [7]. If true, it 
would have opened a big new class of phenomena and 
novel group materials with quite nontrivial properties. 

However there are also arguments [ 8 , 9], that the 
situation with quantum effects in orbitals may be 
not so simple and not exactly analogous to that in 
spin systems. It is predominantly connected with the 


intrinsically strong orbital-lattice coupling, as a result 
of which orbital degrees of freedom become ’’heavy”, 
essentially classical (or one needs to treat also lattice 
vibrations quantum-mechanically, as is done in vibronic 
physics). Specihcally, for YVO 3 Fang and Nagaosa [10] 
proposed an alternative, essentially classical explanation 
of the experimental findings of [7], thus casting serious 
doubts on the importance of orbital quantum effects in 
RVO 3 - practically the only real systems for which these 
effects were claimed to be observed. 

To clarify this important issue, we carried out a 
detailed experimental study of a material of the same 
group and with similar properties as YVO 3 - LUVO 3 , 
using polarized neutron and resonant x-ray scattering. 

Specifically, the intriguing proposal made for the 
RYO 3 compounds was the existence of a novel state - 
an orbital-Peierls dimerization in the c-direction [ 6 , 7]. 
This proposal was based on the observation of a gap 
in the spin-wave spectrum along L, which implied an 
alternation of exchange interaction between consecutive 
layers in the c-direction. The alternative proposal of [10] 
used the more conventional picture of JT distortion and 
orbital order ( 00 ) alternating between consecutive ab- 
planes. However, both these theories produce virtually 
identical spinwave dispersions along the c-direction (see 
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(a) Phase II,T=95 K (c) 



FIG. 1: Magnetic ground state configurations (a,b), spin 
waves (c,d) and relevant exchange coupling schemes for the 
two magnetic phases of LuVOa. The magnetic (orbital) struc¬ 
tures are derived from neutron polarimetry (resonant x-ray 
scattering) data while the spin waves were measured using 
neutron spectroscopy. Calculations shown in (c) are for the 
Jahn-Teller model (solid lines with line width indicating the 
calculated spectral weight) and orbital-Peierls model (dashed 
line), as discussed in the text. Colors indicate the different 
polarisation states. 


Fig. Ic, left part where both theories coincide). The 
spin-wave behaviour of the two competing models can 
be drastically different, but only for carefully selected 
reciprocal lattice directions. Our results along the 
[0K2] direction (Fig. Ic, right part) conclusively favor a 
model based on Jahn-Teller distortions rather than the 
orbital Peierls dimerisation for RVO 3 . Another limiting 
experimental factor in order to distinguish between 
these two pictures is the fact that the material exhibits 
structural, magnetic and orbital transitions at the same 
temperature. Therefore, spin-polarized neutron and 
resonant x-ray beams are essential in order to distinguish 
between magnetic, orbital and structural signals. This 
is of utter importance since a prerequisite to answer 
the main important question is to firstly determine 
the precise magnetic structure. Based on that one can 
make a quantitative analysis on the magnetic-excitation 
intensities of the relevant phase. 

LuVOa is the end-member of the HVO 3 family and 
has the smallest ionic radius, giving rise to increased 
octahedral distortions. We therefore expect it to be 
most sensitive to a possible orbital-Peierls dimerisation. 
The ion {S = 1) is in an octahedral environment 


of O atoms (within the perovskite structure) with the 
two 3(J electrons occupying t 2 g“Orbitals. One electron 
occupies the xy orbital, and the other - one of the 
two doubly-degenerate xz or yz orbitals (or their linear 
superposition). In LuVOa the interplay between spin 
and orbital physics is responsible for the rich phase 
diagram indicated by bulk measurements [ 11 , 12 ], as for 
the well studied YVO 3 [7, 13-19]. Upon cooling, LuVOa 
first enters an orbitally ordered phase at Too = 177 K 
(phase I) followed by magnetic ordering at Tsoi = 105 
K (phase 11). Below that, yet another orbital-magnetic 
phase transition takes place, at Tso 2 = 82 K (phase 
III) [11]. This information, together with the phase 
numbering used throughout this text, is given to the top 
of Fig. 2. 

Experimental setup details can be found in the Sup¬ 
plementary Information. The resonant x-ray scattering 
(RXS) experiment, being sensitive to anisotropic prop¬ 
erties of the tensorial cross-section, yields the charge 
forbidden Bragg reflections arising from the 00 and 
shown in Figs. 2(a), 2(b). Specifically, G-type orbital 
ordering is revealed in phases I+II, while G-type is 
found for phase III (schematic diagrams of the order 
in Fig. 2). Ab-initio calculations show that electric 
dipole transitions dominate the cross-section in this 
case. The neutron data shown in Figs. 2(c) and 2(d) 
complement the x-ray data; they show primarily G-type 
spin ordering for phase II (but with a small admixture 
of G-type) and G-type spin ordering only for phase III. 
This combined neutron-RXS result is in agreement with 
the Goodenough-Kanamori rules, see [20]. 

These measurements demonstrate that LuVOa is an 
antiferromagnet (AF) with fe=0. This, combined with 
structural phase transitions occurring at the magnetic 
ordering temperatures, necessitates the use of polarized 
neutrons to determine the magnetic structures in detail. 
Spherical neutron polarimetry is used to measure polar¬ 
ization matrices for selected reflections in both magnetic 
phases, as discussed in the supplementary information. 
The unpolarized intensity data (E5) were used for 
cross-check as well as for normalising the moment size. 

Phase III has a collinear G-type magnetic structure 
as plotted in Eig. 1(b), consistent with the Pbnm 
orthorhombic space group. The magnetic moments are 
pointing purely along the crystallographic G-axis. Eor 
phase II, beam depolarisation arising due to orientation 
domains implies that the space group can no longer be 
Pbnm. Furthermore, in order to have G-type orbital 
ordering, as observed from our resonant x-ray data, it is 
necessary to lose the mirror plane perpendicular to the 
c-axis, also inconsistent with Pbnm. We could thus fit 
our data with the lower monoclinic P2i/b space group 
and a canted magnetic structure, as shown in Fig 1(a). 

In order to gain deeper insight into structural details 
affecting the precise environment of the ions, 

high-resolution powder x-ray diffraction experiments 
were performed (see suppl. material for details). The 
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FIG. 2: Spin and orbital order parameters of LuVOa as a 
function of temperature, revealed by (a,b) x-ray and (c,d) 
neutron scattering. The complementarity of the techniques 
gives a clear picture of the ordered structures: (a) G-type or¬ 
bital order for phases I, II and (b) C-type orbital order for 
phase III. At the same time, we observe (c) G-type spin order 
for phase III while (d) G-type spin order for phase II, in ac¬ 
cordance with the Goodenough-Kanamori rules [? ]. Note the 
small G-type magnetic component in phase II, which gives an 
overall canted structure. An overview is schematically drawn 
on top of the figure. and refer to spins ”up” and 

"down” or an orbital configuration, for instance "dyz” and 
"dji,” respectively. 


orbitals plotted in Fig. 1 follow precisely the octahedral 
tilts and distortions as determined from refinements of 
these high-resolution x-ray data. 

The spin wave dispersion was measured in both 
phases to help modelling the microscopic Hamiltonian. 
Polarized neutrons were used to measure the dispersion 
in phase II (shown in Fig. Ic), in order to disentangle 
the spin waves from phonon modes, which show similar 
dispersions. Further, we were able to separate magnetic 
excitations with different polarization states and 


standard Blume-Maleev coordinate system nota¬ 
tion [21]), which is important in the subsequent analysis. 
Magnons in phase III were measured by unpolarized 
neutron spectroscopy. 

The spin wave dispersion was modelled using 

linear spin-wave theory (LSWT) and the SpinW 

library [22] with a simple Hamiltonian of the form 

H = E JijSiSj + Ham where Han is the usual easy- 
<i,j> 

axis single-site anisotropy term like —KS^, but with the 
local easy axes different for different sites, see below. 
Jij exchange parameters are shown in Fig. 1 (positive 
Jij denotes AF coupling). Phase HI is fitted with 
two Heisenberg exchange parameters and a single-ion 
anisotropy term (parameters given in Table 1), yielding a 
perfect agreement with our data, as shown in Fig. 1(d). 
The single ion anisotropy term gives rise to two split 
modes as well as the overall gap of the system, found to 
be Egap = 3.7(1) meV. The situation is more complex 
and very different, for the intermediate temperature 
phase H. Magnons are observed at both Q = (0 1 2) 
and Q = (0 1 1) due to the canted structure {C + G 
types). This new phase is gapped, with four modes and 
a splitting of «5 meV between the two sets of branches. 
The c axis energy scale is reduced by more than a factor 
of two while the doubling of the number of modes is in 
accordance with two inequivalent V sites (lowering of 
the space group). 

In order to determine the Hamiltonian, in this more 
complex case, it is necessary to model data along a 
further direction (& axis. Fig. 1(c)). The dispersion along 
the b axis rises to 27 meV, an energy which can only be 
explained by assuming the interaction scheme of Fig. 
1 (a) with two distinct ab plane couplings alternating 
along c, as well as a ferromagnetic coupling along the c 
axis [10]. The parameters of this Hamiltonian which fit 
the data very precisely are also given in Table 1. In order 
to account for the canting of the moments, two individ¬ 
ual directions were assumed for the single ion anisotropy 
term: easy b and c axes for Vi and V 2 respectively, 
located in alternating a6-planes in accordance with 
the P2i/b space group symmetry and the interaction 
scheme. Note that this model accounts not only for the 
measured spin wave energies and intensities, but also 
correctly indexes the polarization states and 


TABLE I: Exchange parameters for both magnetic phases of 
LuVOa. Positive sign corresponds to AF coupling. 


Phase III (meV) 

Phase II (meV) 

Jab = 4.24(21) 

Jc = 5.95(19) 

K = (-0.48(12), -0.06(2), 0) 

Jabi = 0.82(3) 

Jab2 = 5.99(3) 

Jc = -1.29(2) 

Ki = (0, 0.66(5), 0) 
K 2 = (0, 0, 0.66(5)) 
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FIG. 3: Magnetic scattering intensity as measured in the in¬ 
termediate phase of LuVOa {T—95 K). The directions shown 
are the c* and 6*-axes (left and right panels respectively). 
The exchange models being compared are either the 3D Jahn- 
Teller model (solid red lines) or the orbital Peierls model (dot¬ 
ted black lines). 


further respects the magnetic structure as well as the 
symmetry of the P2ijb space group. An alternative 
orbital Peierls model with alternating exchange along 
the c-direction and a single coupling in the ab plane 
would require a lowering of the space group due to 
loss of a mirror plane along the crystallographic c axis. 
This model [7] was fitted to the energies of our dataset, 
but predicts an almost flat mode at 15 meV for the 
dispersion along the b axis (see black dashed line in Fig. 
Ic) with zero intensity, contrary to our observation. The 
direct comparison of both models with the data (Fig. 
3) shows that indeed the coupling scheme of Fig. 1(a) 
(continuous red lines) is in excellent agreement with 
our data. In contrast, the orbital-Peierls model with 
exchange couplings alternating along c (dashed black 
lines) shows large overall discrepancies in the spectral 
weight distributions. The figure shows the calculated 
intensities after fitting the dispersion energies (Fig. 1), 
with a single scaling factor. Note that the symmetry 
of the monoclinic P 2 i /6 space group is consistent with 
inequivalence of the bonds in alternate ab planes (3D 
JT model), but incompatible with inequivalence of 
the c axis bonds between these planes (orbital-Peierls 
scenario). This is in addition a very strong symmetry 
argument in favour of the former model. 

Comparing phase III of LuVOa at low temperature 
with YVO 3 we find a direct agreement between the 
ground state and model Hamiltonian as discussed in 
Ref. [7]. The moment direction, inelastic gap-size and 
band-width match well. The main difference is a 30% 
smaller ordered moment and variations of up to 40% in 
the exchange couplings and anisotropy parameters which 
can be attributed to the details of lattice distortions 
and exchange pathways in each case. However, the 
experimentally determined magnetic structure of phase 


II at intermediate temperature differs from all other 
RYO 3 members [7, 13, 23, 24]. The use of polarized 
neutron diffraction enables precise determination of the 
moment direction (fee-plane), due to the k = 0 magnetic 
propagation vector combined with structural distortions 
across the magnetic phase transitions. Furthermore, by 
using full polarization analysis it is possible to predict 
the subtle lowering of the space group from orthorhom¬ 
bic to monoclinic, based on symmetry arguments and 
observation of orientation domains. 

By employing polarized neutron spectroscopy and, in 
particular, by making extended measurements in two 
inequivalent directions in reciprocal space, we were able 
to distinguish between model Hamiltonians in which 
either the exchange Jat or Jc alternate along the c-axis 
between two values (see Figures la, c and 3). The 
clear evidence for the first model {Jab alternation) is 
the superiority of dispersion and intensity fits to this 
model, compared to the other (J^ alternation). On 
top, this is in excellent agreement with ab — initio 
theoretical calculations given in Ref. [10]. The alternate 
Jab exchange parameters calculated for YVO 3 are 0.8 
and 5.3 meV, and our fitted results for LUVO 3 are 0.8 
and 6 meV, strongly supporting our conclusion. These 
first-principle calculations [ 10 ] are based on the precise 
experimentally determined JT distortions of YVO 3 and 
related vanadates [14-16, 25] which do indeed alternate 
between adjacent afe-planes, exactly in the same fashion 
as Jab in our Hamiltonian. The Jab exchange parameter 
is very sensitive to these JT distortions and orbital 
ordering, because it depends on subtle competition 
between various exchange processes [10]. This explains 
the large difference between the two alternating values 

of Jah- 

Based on these new findings, we conclude that the or¬ 
bital fluctuations -which are inherent to these systems- 
are in fact suppressed by the JT distortions. In LuVOs 
this results in an overall 3D spin-orbital structure, rather 
than a quasi-lD orbital dimerized chain. 

Summarizing, we have carried out a detailed analysis 
of the interplay amongst spin, orbital and lattice degrees 
of freedom in the Mott insulator LUVO 3 . By combining a 
variety of experimental methods, we are able to uniquely 
determine the magnetic and orbital states and to model 
the spin Hamiltonian in the two magnetic phases. These 
results show that the features, attributed previously 
to an orbital-Peierls state (’’orbital-singlet”, similar to 
spin-singlet dimers) which could have appeared due to 
quantum effects in orbitals, are in fact a consequence of 
the static orbital ordering and corresponding Jahn-Teller 
distortion. Yet, we cannot rule out that the orbital 
quantum fluctuations may still be present in some form, 
maybe in different materials. This question deserves 
further study. 
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The measurements were performed on both powder 
samples and single-crystal fragments using synchrotron 
x-rays to determine the crystal structure and orbital or¬ 
der pattern. The instruments used were the X04SA-MS 
beamline at SLS, Villigen [1] the 6ID-B resonant x-ray 
beamline, APS, Argonne. The diffraction part was per¬ 
formed using A = 0.620 A. For the resonant x-ray exper¬ 
iment the vanadium Al-edge {E = 5.482 keV) was used 
in order to acess the orbital order wavevector transfer 
Q required. In order to probe the magnetic order pa¬ 
rameter and excitations directly we have used diffraction 
and spectroscopy neutron scattering techniques respec¬ 
tively. The four-circle diffractometer E5 at HZB, Berlin 
and the D3 diffractometer at ILL with the spherical neu¬ 
tron polarimetry option (CRYOPAD) were operated at 
A = 2.38 A and A = 0.825 A with a ®He spin analyser 
respectively. The FLEX triple-axis spectrometer (TAS) 
at HZB [2, 3] was operated with a fixed final wavevector 
kf = 1.55 A“®, and a Be-filter to reduce second-order 
contamination. The thermal-neutron TAS INS and IN20 
at ILL, Grenoble (the last with XY Z-polarization analy¬ 
sis) were used with fixed final wavevector kf = 2.662 A“® 
and a pyrolytic graphite filter for second-order contami¬ 
nation. The large single crystal used in the spectroscopic 
measurements weighs 4 g. 

The magnetic structure of LuVOs was determined in 
both phases II and III by using spherical neutron po¬ 
larimetry. The resulting polarisation matrices were mod¬ 
elled with the code Mu fit [4] based on the Blume-Maleev 
equations [5, 6]. Each element of the matrix was averaged 
over the different domains, the populations of which were 
treated as free parameters and fitted by the code. Table 
SI gives the data and the fits for a selection of reflections. 

Phase III is a collinear G—type AF with moments 


TABLE SI: Representative measured and calculated polariza¬ 
tion matrices used to determine the magnetic structures for 
both magnetic phases of LuVOa. 


phase III (low T phase) 


h 

k 

1 

initial P 

Observed 

Calculated 

1 

0 

1 

1 

0 

0 

- 0 . 65 ( 2 ) 

0 . 07 ( 1 ) 

0 . 05 ( 1 ) 

- 0.62 

0.03 

0.00 

1 

0 

1 

0 

1 

0 

- 0 . 00 ( 1 ) 

0 . 98 ( 2 ) 

- 0 . 02 ( 1 ) 

0.00 

1.00 

0.00 

1 

0 

1 

0 

0 

1 

- 0 . 03 ( 1 ) 

0 . 05 ( 1 ) 

- 0 . 62 ( 2 ) 

0.00 

0.03 

- 0.62 

3 

0 

1 

1 

0 

0 

0 . 68 ( 6 ) 

0 ( 1 ) 

0 ( 1 ) 

0.59 

- 0.03 

0.00 

3 

0 

1 

0 

1 

0 

0 ( 1 ) 

0 . 98 ( 8 ) 

0 ( 1 ) 

0.00 

1.00 

0.00 

3 

0 

1 

0 

0 

1 

0 ( 1 ) 

0 ( 1 ) 

0 . 79 ( 7 ) 

0.00 

- 0.03 

0.59 

1 

0 

3 

1 

0 

0 

0 . 99 ( 1 ) 

0 ( 1 ) 

0 ( 1 ) 

0.99 

0.00 

0.00 

1 

0 

3 

0 

1 

0 

0 ( 1 ) 

1 . 01 ( 1 ) 

0 ( 1 ) 

0.00 

1.00 

0.00 

1 

0 

3 

0 

0 

1 

0 ( 1 ) 

0 ( 1 ) 

0 . 99 ( 1 ) 

0.00 

0.00 

0.99 

1 

0 

1 

-1 

0 

0 

- 0 . 63 ( 2 ) 

- 0 . 04 ( 1 ) 

0 . 02 ( 1 ) 

- 0.62 

0.03 

0.00 

1 

0 

1 

0 

-1 

0 

0 . 01 ( 1 ) 

0 . 97 ( 2 ) 

- 0 . 00 ( 1 ) 

0.00 

1.00 

0.00 

1 

0 

1 

0 

0 

-1 

- 0 . 03 ( 1 ) 

0 ( 1 ) 

- 0 . 62 ( 2 ) 

0.00 

0.03 

- 0.62 


phase II (intermediate T phase) 


H 

K 

L 

initial P 

Observed 

Calculated 

0 

1 

0 

1 

0 

0 

- 1 . 00 ( 10 ) 

0 . 05 ( 7 ) 

0 . 01 ( 7 ) 

- 1.00 

0.00 

0.00 

0 

1 

0 

0 

1 

0 

0 . 01 ( 7 ) 

1 . 02 ( 10 ) 

- 0 . 01 ( 7 ) 

0.00 

1.00 

0.00 

0 

1 

0 

0 

0 

1 

- 0 . 02 ( 7 ) 

- 0 . 01 ( 7 ) 

- 1 . 03 ( 11 ) 

0.00 

0.00 

- 1.00 

1 

0 

0 

1 

0 

0 

- 1 . 00 ( 7 ) 

0 . 01 ( 5 ) 

0 . 02 ( 5 ) 

- 1.00 

0.00 

0.00 

1 

0 

0 

0 

1 

0 

- 0 . 00 ( 4 ) 

0 . 11 ( 4 ) 

0 . 06 ( 4 ) 

0.00 

0.10 

0.05 

1 

0 

0 

0 

0 

1 

- 0 . 03 ( 4 ) 

0 . 05 ( 4 ) 

- 0 . 10 ( 4 ) 

0.00 

0.05 

- 0.10 

1 

0 

-2 

1 

0 

0 

- 1 . 00 ( 17 ) 

0 ( 1 ) 

0 ( 1 ) 

- 1.00 

0.00 

0.00 

1 

0 

-2 

0 

1 

0 

0 ( 1 ) 

- 0 . 42 ( 13 ) 

0 ( 1 ) 

0.00 

- 0.41 

0.05 

1 

0 

-2 

0 

0 

1 

0 ( 1 ) 

0 ( 1 ) 

0 . 38 ( 12 ) 

0.00 

0.05 

0.41 

1 

0 

2 

1 

0 

0 

- 1 . 00 ( 15 ) 

0 ( 1 ) 

0 ( 1 ) 

- 1.00 

0.00 

0.00 

1 

0 

2 

0 

1 

0 

0 ( 1 ) 

- 0 . 44 ( 11 ) 

0 . 11 ( 9 ) 

0.00 

- 0.41 

0.05 

1 

0 

2 

0 

0 

1 

0 ( 1 ) 

0 . 13 ( 10 ) 

0 . 37 ( 11 ) 

0.00 

0.05 

0.41 


pointing along the crystallographic c—axis as shown in 
Fig. 1(b) of the main text. This is consistent with the 
Pbnm orthorhombic space group. The moment vector 
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FIG. SI: Variation in the octahedral environment of V®'*' in 
LuVOa as a function of temperature: (a) lattice parameters 
and (b) V-0 bond lengths refined from high-resolution syn¬ 
chrotron x-ray data. 

on each V site is = [0 0 1.20(2)] fis- The rela¬ 

tive populations of the 180° domains that break time- 
reversal symmetry is 48(3)%-52(3)%. Phase II is more 
complex, partly because it combines C— and G—type 
ordering. We find Mq = [0 0.66(5) 0.74(5)] and 
Mq = [0 0.32(8) ±0.36(8)] for the two contributions, 
i.e. the moments are canted and lie in the 6 c—plane. The 
’±’ ambiguity arises from the existence of equally popu¬ 
lated orientation domains, making it impossible to cor¬ 
relate the orientation domains of the G and C modes. 
However, the ’±’ solution results in a moment size larger 
than possible for a free ion. As shown in Fig. 1(a), 
the solution gives a lattice with a uniform overall mo¬ 
ment size per site = 1-10 fiB- Neither the spin 

orientation derived from the polarized neutron data nor 
the combination of G and G modes is compatible with 
space group Pbnm, but rather point to a lowering of the 
crystal symmetry to the monoclinic space group P2i/b. 
This is consistent with the presence of orientation do¬ 
mains (htted populations 47(3)%-53(3)%). The neutron 


polarimetry data and P2i/b space group of this phase are 
consistent with the symmetry of the Hamiltonian mod¬ 
elling our spin-wave data, as discussed in the main text. 

X-ray powder diffraction patterns were collected at 60 
different temperatures and htted using the FullProf pro¬ 
gram [7]. The phase transitions are visible in the unit 
cell parameters plotted in Fig. SI(a). It is most instruc¬ 
tive, however, to inspect the V-0 distances shown in Fig. 
SI(b), which are those dominating the super-exchange in¬ 
teractions. The V-01 (apical) distance stays relatively 
constant as a function of temperature, suggesting that 
the occupancy of the dxy orbital is unchanged. By con¬ 
trast the lengths of the two distinct V-02 bonds change 
quite drastically across the HI—^H phase transition and 
the difference between them changes sign pointing to a 
corresponding swap in the occupancies of the dxz and dyz 
states. Their absolute difference is minimised in phase H 
and changes sign at T = 100 K. We can therefore ex¬ 
pect drastic changes in the effective Hamiltonian across 
the transition and orbital fluctuations in phase H, as pre¬ 
dicted by theory [ 8 ]. Structural refinements of our high- 
resolution x-ray data with the two space groups {Pbnm 
or P2i/6) give indistinguishable fit quality and x^. In 
other words, x-ray diffraction alone is unable in this par¬ 
ticular case to solve the crystallographic structure, even 
with the highest possible resolution available since the 
monoclinic splitting is too small. Our case highlights 
the power of polarized neutron data in combination with 
symmetry analysis. 
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